Inelastic scattering of atoms in a double well 
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We study a mixture of two light spin- 1/2 fermionic atoms and two heavy atoms in a double well 
potential. Inelastic scattering processes between both atomic species excite the heavy atoms and 
renormalize the tunneling rate and the interaction of the light atoms (polaron effect). The effective 
interaction of the light atoms changes its sign and becomes attractive for strong inelastic scattering. 
This is accompanied by a crossing of the energy levels from singly occupied sites at weak inelastic 
scattering to a doubly occupied and an empty site for stronger inelastic scattering. We are able to 
identify the polaron effect and the level crossing in the quantum dynamics. 



I. INTRODUCTION 

Ultracold bosonic and fermionic atomic gases in optical lattices can be used as a toolkit for the investi- 
gation of fundamental condensed matter physics models [1] . Recent experimental work opened a field to 
study quantum states in optical lattices, such as superfluid and Mott states [2, 3], where the interparticle 
interaction can be controlled by a magnetic field via a Feshbach resonance [4]. Spin-dependent effects 
[5-9], frustrated spin systems [10], the formation of dimers from fermionic atoms [11, 12], and mixtures of 
two atomic species [13-15] provide opportunities for creating and studying even more complex quantum 
states. 

The optical lattices are robust and free of phonons. On the other hand, the electron-phonon interaction 
in a solid leads to a rich physics. It is important for superconductivity, the Peierls instability, polaron 
effects and many other phenomena. With more progress in the atomic and laser physics, the coupling of 
ultracold atoms in an optical lattice to bosonic degrees of freedom may be achieved and thus can mimic 
the dynamics of electrons in the presence of phonons. Recently, ultracold atoms confined to an optical 
resonator were proposed to study the effect of coupling between the atoms and the photons field, which 
leads to an effective Hubbard Hamiltonian with long-range interaction [16] and to an interesting phase 
diagram [17]. Bosc- Fermi mixtures can also provide an insight into the role of bosons in the dynamics of 
fcrmions, where the condensed bosons lead to fermionic pairing [18] and charge density waves of fcrmions 
[19]. An interesting example of the latter are dimer states. They have been discussed in solid-state 
systems [20], in the Holstein-Hubbard model [21] and recently also for an ultracold Bose gas with ring 
exchange [22]. 

More recently, ultracold gases were employed to study the dynamics of quantum states, including the 
"collapse and revival" behavior [23]. Here it is important to distinguish between small systems with 
a few atoms and many-body systems with a large number of atoms [1]. For instance, it was observed 
experimentally that in a small system with two spin- 1/2 atoms the spin dynamics and the particle 
dynamics are completely separated, similar to the spin-charge separation in one-dimensional systems 
[24, 25]. Another example for a restricted dynamics in small atomic systems are entangled squeezed 
states in a Bose-Einstein condensate, whose atoms are distributed over a small number of lattice sites 
[26] . Both observations indicate that the dynamics of small atomic systems can be restricted to a subspace 
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of the entire Hilbert space available for the model Hamiltonian. This can also mean that the system never 
reaches the groundstate of the Hamiltonian if it was prepared in an excited state, reflecting the fact that 
the initial state has no overlap with the groundstate. For two spin-1/2 atoms in a double well, described 
by a Hubbard model, this is a direct consequence of the fact that the eigenstates do not mix pairs of singly 
occupied sites with pairs of empty and doubly occupied sites [27]. On the other hand, mixing of these 
states in a macroscopic system, enforced by inelastic scattering with other atoms, can lead to a first-order 
quantum phase transition from singly occupied sites to doubly occupied sites. This was observed in a 
mixture of light and heavy atoms, where the latter are in a Mott state [21]. This case can be described 
by a Bose- Fermi model that is known in solid-state physics as the Holstein-Hubbard (HH) model [28]. 
Adjusting physical parameters, such as the optical-lattice parameters (frequency and amplitude of the 
Laser field) and the fcrmion-fcrmion interaction through a Fcshbach resonance, enables us to prepare such 
a system not only in the ground state but also in its excited states and to study its dynamics. Although in 
a small system there is no phase transition for the ground state, dynamical properties of excited states can 
change qualitatively due to inelastic scattering described by the HH model. Among other effects, there 
is renormalization of the tunneling rate, known as the polaron effect, which was recently also discussed 
for an ultracold Bose gas [29] . 

Our interest is to study the effect of inelastic scattering of two spin-1/2 fermionic atoms in a double- 
well potential with repulsive Hubbard interaction and an additional scattering with a heavy atom in each 
potential well. Each of the heavy atoms have a harmonic oscillator spectrum and can exchange energy 
with the tunneling spin-1/2 atoms. 

In this paper, we study a double well potential filled with heavy (bosonic or fermionic) atoms (HA) (e.g. 
87 Rb or 40 K), one at each well (see Fig. 1). A double-well potential can be realized by superimposing 
two periodic potentials with different periodicity and the dynamics of one or two atoms can be easily 
studied. The tunneling of these HA is neglected since the potential barrier between the wells is sufficiently 




FIG. 1: Two LFA in a double well potential with one HA in each well. The HA are harmonic oscillators and are 
represented by horizontal lines (energy levels). The LFA can tunnel from one well to the other one. The wells 
can be singly occupied (left panel) or one well can be doubly occupied and the other one empty (right panel). 

high. Excitations are only due to collisions with other atoms. For this purpose two light fermionic atoms 
(LFA) (e.g. 6 Li), prepared in two hyperfine states denoted as | I), are added to the system. These 
atoms can tunnel because of their low mass and scatter by the HA. It is assumed that the HA experience 
a harmonic potential of a well, at least at low energies. Then their excitations are harmonic-oscillator 
states. During the scattering process the HA can also transfer energy to the LFA. Moreover, the light 
fermions experience local (on-site) repulsion. 

The paper is organized as follows. The Holstein-Hubbard model is introduced and discussed in Sect. 
II. In Sect. Ill we introduce a restricted model with at most one phonon excitation per well. Then in 
Sect. IV the effective Hamiltonian of the unrestricted Holstein-Hubbard model is defined, its spectral 
properties are studied and compared with those of the restricted model of Sect. III. Based on this 
effective Hamiltonian we study the dynamics of the quantum states in a double well, the spectral density 
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and the spin imbalance in Sect. V. 



II. HOLSTEIN-HUBBARD MODEL 



The atomic mixture of LFA and HA can be well described by the Holstcin-Hubbard model [28] : 

H = -J ]T ]T cj ffCj , CT + h.c. + J2 [wo&]&j + 9(b] + &j)(njt + «jj + Un n n u ] . (1) 

The first term describes the tunneling of LFA with spin a (=^, I) between nearest-neighbor wells. These 
are defined by fermionic creation and annihilation operators c. and c- }ai respectively. The HA form a 
Mott state and are presented as harmonic oscillators at each well with eigenfrequency u>o, assuming that a 
HA in one well is excited independently of the HA in the other well. Thus they can be considered as local 
phonons and are described by the bosonic creation and annihilation operators bj and bj. The phonons 

couple to the light atoms with strength g ~ (elT^I/), where V is the interaction between LFA and HA, 
j/) denotes the ground state of a HA, while |e) denotes its first excited state. The fourth term describes 
the interaction between two LFA at the same well, where U is a local repulsive interaction between the 
LFA 

This lattice model describes the quantum phase transition in a half-filled system from singly occupied 
lattice wells (Neel state) to a mixture of doubly occupied and empty wells (dimer state). Now we restrict 
the lattice model to the two sites of the double well potential, choosing the coordinates j = 1, 2 for the 
wells. Ignoring the tunneling of the LFA and applying a unitary transformation to the remaining part 
of the Hamiltonian we can decouple fermionic and bosonic degrees of freedom and get the transformed 
local Hamiltonian [21] 



j = l,2 



g2 



(2) 



with 7 = U — 2g 2 /ojq is the effective Hubbard coupling. For a system with two fermions the ground state 
energy of H~ ( is given by 

E = e - 73 , (3) 

with e = U — Ag 2 /loq and 2j is the number of singly occupied wells. The coupling 7 controls two different 
regimes: for 7 > 0, the ground state has two singly occupied wells (j = 1) and energy —2g 2 /uiQ, while for 
7 < 0, there are one doubly occupied well and one empty well (j = 0) and the energy is e. Thus at 7 = 
there is a transition, where the system changes from two singly occupied wells to a doubly occupied well 
and an empty well. Moreover, both states are degenerate even for 7 7^ because the local Hamiltonian 
Hj does not determine how the spins, and the empty and doubly occupied wells are distributed in the 
double- well potential. Tunneling in the Hamiltonian H lifts these degeneracies [28]. 



III. DOUBLE WELL WITH AT MOST ONE PHONON EXCITATION PER WELL 

Even for a double well with two fermions the Hilbert space of the Hamiltonian in Eq. (1) has infinite 
dimensions due to the phonon excitations. For nonzero tunneling rate J this becomes a difficult problem. 
However, to study quantitatively the effect of inelastic scattering we can restrict the phonon excitations. 
The simplest case is a model in which each of the localized atoms is a two-level system, similar to the 
Jaynes-Cummings model in quantum optics [30]. For J = the four lowest eigenvalues are degenerated 
in pairs, as shown in the left panel of Fig. 2. The ground state has a cusp, but the latter does not coincide 
with the exact ground state for J = given in Eq. (3). A nonzero tunneling rate J lifts the degeneracy 
and all four eigenvalues are distinct now (cf. left panel of Fig. 3). We notice that for large values of g 
the curves do not merge (i.e., the distances between the first and the third, and the second and fourth 
curves, respectively, remain nonzero). 
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FIG. 2: Model without tunneling (J = 0): Four degenerate eigenvalues of the two-sites problem with at most one 
phonon excitation (left panel) and the four lowest poles of the resolvent of Eq. (4) with the effective Hamiltonian 
H e ff in Eq. (5) (right panel) for (7 = 2 and ojo — 50. Curve 1 (curve 2) represent levels with two singly (one 
empty and one doubly) occupied wells. Left panel: The two curves cross at g ~ 7.6. The energy of the ground 
state has a cusp, but it does not coincide with the exact ground state for J = given in Eq. (3). Right panel: 
The two curves cross at the smaller value g = 5\/2 « 7.07. 




FIG. 3: Model with tunneling (J = 1): Four eigenvalues of the two-sites problem with at most one phonon 
excitation (left panel) and the four lowest poles of the resolvent of Eq. (4) with the effective Hamiltonian H e g 
of Eq. (5) (right panel) for U — 2 and ojo = 50. Left panel: The curve 2 and the curve 3 cross at g ss 7.6. For 
very large values of coupling g the distances between curve 1 and curve 3, as well as between curve 2 and curve 
4, remain nonzero. Right panel: The curve 2 and the curve 3 cross at g m 7.1. For very large values of coupling 
g the distances between curve 1 and curve 3, as well as between curve 2 and curve 4, vanish as ~ r ' 2 /\j\- This is 

_ 2 / 2 

due to the polaron effect, which accounts for rescaling of the single fermion tunneling J — > r = Je 9 /UJ ° . 



IV. EFFECTIVE HAMILTONIAN FOR MANY PHONON EXCITATIONS 



Now we consider the full Holstein- Hubbard Hamiltonian in Eq. (1) and treat its spectral properties in an 
approximative manner. The main idea is to study the evolution of the quantum system \^t) = e~ ltH \^o), 
beginning with the initial state |Wo). The evolution is a walk through the entire Hilbert space that is 
accessible for the Hamiltonian. The Recursive Projective method (RPM) organizes this walk by projecting 
iteratively on a sequence of subspaces. The main advantage of this method is that the walk visits each 
subspace only once [27, 31, 32]. The approximation method within the RPM consists of ignoring some 
part of the Hilbert space that contributes with a low probability to the walk and leads to an effective 
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Hamiltonian. Details of the application of the RPM to the Holstein-Hubbard Hamiltonian can be found 
in Ref. [21]. In the following we start from the effective Hamiltonian that was derived in Ref. [21], 
to study the dynamics of the LFA in the double-well potential. The advantage of this method is that 
it enables us to study the effect of finite tunneling of the LFA as well as arbitrary number of phonon 
excitations. 

In order to derive an effective Hamiltonian, we project the full Hilbcrt space of the Hamiltonian in 
Eq. (1) onto the Hilbert space spanned by the four Fock states \*t,i) , I > |0,TI) an d |tl;0) with a 
projector Pq, such that the resolvent of the projected Hamiltonian is given by 

G (z)=P (z~H)- 1 P = (z-H^z))' 1 . (4) 

The effective Hamiltonian for the double well can be evaluated recursively by the RPM. Under the 
assumption that the tunneling rate of LFA J is small compared to the other parameters of the system 
(e.g. J <C ojq, J -C £/ and J <C g) the recursion relation can be truncated, which gives [21] 

H e s(z) RJ ^ [- Tc lcr C 2<? + K l( z )c\<T C 2° C 2cr> C lc' + K 2 ( Z ) C \a C ^ C la' + hx - + \ (. n lv n lcr' + n%„n% a i ) 

(5) 

where the indices 1 and 2 represent the left and right sites of the double well, respectively. This Hamil- 
tonian describes three different tunneling processes, namely the tunneling of single fermions with rate r 
(first term), the exchange of spins with rate K\(z) (second term), the tunneling of fermionic pairs with 
rate K 2 (z) (third term) and the on-site interaction between fermions with strength 7 (fourth term). The 
tunneling rate J of single fermions in Eq. (1) is now renormalized as 

r = e -9 2 /^oJ, (6) 

which is the well-known polaron effect [28]. The spin-exchange parameter K\ and the pair tunneling 
parameter K 2 are given by the expressions [21] 



ml z — 2e + 27 — worn 



and 



L — ' m! z — le — LUnm 

m—l 

In order to avoid the singularities of the coefficients K\ and K2, we assume that loq 3> U, g, which is valid 
for a deep and tight double-well potential. 

The energy levels of the system is given by the poles of Eq. (4). Thus the variable z is fixed by solving 
the equation 

det[z - fleff(*)] = 0. (9) 

To solve this equation we first diagonalizc the 4x4 effective Hamiltonian H e ff for a fixed parameter z 
and find its eigenvalues Xj(z) (j = 1,...,4). Then we solve z = Xj(z) for each of the four eigenvalues 
Xj(z) to determine the poles of the resolvent in Eq. (4). An eigenstate \E) (with H e g\E) = E\E)) can 
be written as a linear combination of the four Fock states as 

\E) =a 1 (z)^,i} + a 2 {z)\U) + a 3 (z)\0,n} + a i {z)\n,0) , (10) 

where the coefficients dj(z) run over all possible poles z. In this Fock-state basis the Hamiltonian in Eq. 
(5) reads 



K\ + 6 — 7 —K\ —t —t 

T T 

vi + e K 2 
K 2 K t +e 



H eS (z)=\ T T Kl+e K2 (11) 



G 



whose eigenvalues and the coefficients of the corresponding (non-normalized) eigenvectors in the form of 
Eq. (10) are for doubly occupied lattice sites 

Ai (z) =K 1 {z)-K 2 {z) + e: ai =a 2 = 0; a 3 = -1 ; a 4 = 1 (12) 

and for singly occupied lattice sites 

A 2 (z) = -7 + e: a\ = a? = 1 ; a 3 = a 4 = . (13) 

There are also states with a mixture of singly and doubly occupied sites: Using S = 
^16-r 2 + (-K 1 (z) + K 2 {z) + 7) 2 we have 

X 3 (z) = i (3Kj.(z) + K 2 (z) - 7 - S + 2e) (14) 



with 



and 



with 



-A'i(z) + /r 2 (z)+ 7 + 5 

«i = ; = -a 2 \ a 3 = a 4 = 1 (15) 

4r 



A 4 (z) = 1 (3^x(z) + X 2 (z) - 7 + 3 + 2e) (16) 



-A"i(z) + K 2 (z) + 7 -S 

Oi = = -a 2 ; a 3 = a 4 = 1 . (17) 

4t 

In the right panels of Figs. 2 and 3 we plot the four lowest poles of the resolvent of Eq. (4) with the 
effective Hamiltonian H e s of Eq. (5) as functions of g. In particular, in Fig. 3 curve 1 represents a 
solution of z = X 3 (z), curve 2 a solution z = A 2 (z), curve 3 a solution of z = Ai(z), and curve 4 a solution 
of z = A 4 (z). These four poles are compared with the four eigenvalues of the restricted model with at 
most one phonon excitation of Sect. Ill, shown in the left panels of Figs. 2 and 3. 

In case there is no tunneling (i.e. J = 0) we get degenerate states and only two different eigenvalues 
which cross each other are available (sec Fig. 2). The ground state thus has a cusp at g = ps 7.1 

and coincides with the exact ground state of the Holstein-Hubbard model for vanishing tunneling given 
in Eq. (3). This was not the case, when we considered the exact solution with only one phonon excitation 
in the previous section. 

A nonzero tunneling lifts the degeneracies and leads to a unique ground state. The eigenvalues for 
nonzero tunneling and for U = 2, loq = 50 are shown in Fig. 3. The two lowest excited states still cross 
at around g w 7.1, while the ground state is unique and thus the system does not exhibit the transition 
discussed in Sect. II. As a consequence of the polaron effect, the renormalized tunneling rates r, Ki, 
and K 2 vanish for large g. This implies for the eigenvalues the asymptotic behavior 

Aj ~ A 4 ~ -4g 2 /uj , A 2 ~ A 3 - -2g 2 /oj . (18) 

This is also visible on the right panel of Fig. 3, while for the exact two sites problem with at most one 
phonon excitation (previous section), the eigenvalues do not merge for large coupling g (cf. left panel of 
Fig. 3). 

We plot the coefficients dj (see Eq. (10)) for the ground state of the effective Hamiltonian in Fig. 4. 
There is a crossover from the domination by the singlet states | f, |) (at small coupling g when the effective 
interaction is repulsive, 7 > 0) to the domination of the doubly occupied states (at larger coupling g when 
the effective interaction is attractive, 7 < 0). 

To discuss the consequences of the crossing on observable dynamical quantities we study in the next 
section the dynamics of the quantum states and the spin imbalance of the LFA. This investigation includes 
the spectral density of the model. 



r 




FIG. 4: Coefficients (see Eq. (10)) for the ground state of the effective Hamiltonian H c ff in Eq. (5) (i.e. curve 
1 in the right panel of Fig. 3) with J = 1, U = 2 and cjo = 50. They cross very close to the crossing point 
of the eigenvalues in Fig. 3. The domination of the singlet states 1 1, -l) at small couplings g is changed by the 
domination of the dimer states for large values of g. 



V. DYNAMICS IN A DOUBLE- WELL POTENTIAL 



The description of the dynamics of our quantum system is based on the knowledge of the energy levels, 
the initial quantum state and the overlap of the energy eigenstates with the initial quantum state. In 
other words, if the system with energy levels Ej is prepared initially in state l^o)- its quantum state \^t) 
evolves in time as 

=e- itH \* ) =J2e- UEi \Ej)(Ej\*o) ■ (19) 

3 

The energy levels Ej and the overlap with the initial state (Ej\^o) can be described by spectral density. 
This will be discussed in the next section. 



A. Spectral density 



The return probability of the system to the initial state is calculated from the inverse Laplace transform 
of the projected resolvent of Eq. (4) as [27] 

/dz 
e-^^^z-H^z))^,)— , (20) 

where the contour T includes all the poles of the resolvent. The many-body spectral density is then given 
by 

ps(E) = --Im(y \(z-H cS (z))o 1 \yo) , (21) 

7T u 

where z = E + iS and 5 <^ 1 . For a finite dimensional Hilbcrt space it is a rational function with poles 
3 = ^(7 = 0,1,2,...) 

Ps(E) = 5^ {E _ E ^ + 62 . (22) 

This expression represents Lorcntzian peaks at positions Ej whose heights are \(^>o\Ej)\ 2 /S. Plotting 
ps{E) as a function of E, we can clearly identify the poles Ej of the resolvent Go (see Eq. (4)) and the 
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overlap between the energy state \Ej) and the initial state l^o)- Here we calculate the spectral density 
for the initial state l^o) = I til) f° r different coupling g. Then it should be noticed that the initial state 
is singly occupied and has no overlap with the doubly occupied eigenstate of Eq. (12). The results are 
shown in Fig. 5 for low energies. We observe three peaks, which correspond to the energies shown in Fig. 
3. The central peak represents the dominant energy level for the dynamics. We notice the absence of one 
peak (corresponding to curve 3 in Fig. 3), since the state (|tl>0) + |0,tl))/v2 can not be reached from 
the initial state = I til)- We also observe that before (g = 4) and after (g = 10) the crossing of the 
eigenvalues, there is one dominant central peak and two lower peaks. The characteristic frequencies of the 
dynamics are the differences of the energies of the central peak and the other ones. Close to the crossing 
point (g = 7.06), the two other peaks are symmetric with respect to the central one. Consequently, only 
one frequency appears in the dynamics at the crossing point. This will also be seen in the spin imbalance 
of the next section. For big values of g we observe two large peaks which are very close and a small peak 
far from the central peak. This implies a dominating small single frequency, as also found in the spin 
imbalance (cf. Fig. 8). 
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FIG. 5: The spectral density in Eq. (22) for the initial state |*o) = 1 1- 1) and for J = 1, U = 2, lu = 50, S = 0.02 
and different values of the coupling g. We observe three peaks because the fourth eigenstate is orthogonal to the 
initial state. The central peak represents the dominant energy level and for big values of g the contribution of 
one peak is very small. The frequencies of the spin imbalance are given by the difference between the dominant 
energy level and the other ones. 



B. Spin imbalance 

A recent experimental study of the dynamics of two spin- 1/2 atoms with strong repulsion in a double 
well has revealed that, using two singly occupied wells as the initial state, the single occupation is static 
while the spin oscillates periodically between the two wells with two characteristic frequencies [24, 25]. 
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FIG. 6: Spin imbalance Ni^(t) without inelastic scattering (g = 0, left panel) and with strong inelastic scattering 
(g = 10, right panel) for J = 1, U — 2 and ojo — 50. There are two modes with frequencies uii « 3.24 and 
u)2 ~ 1.24 (left panel) and lui ~ 3.17 and ui2 ~ 1-16 (right panel) contribute to the dynamics. 



This observation has been interpreted by an effective dynamics based on the Heisenberg model. The latter 
can be understood either within a strong-coupling approximation of the underlying Hubbard model [25] 
or within the recursive projection method for a general coupling [27]. Experimentally this has been seen 
by measuring the spin imbalance between the two wells 



In our Holstein-Hubbard model we can vary the coupling g between the LFA and the HA to realize an 
additional interaction. We have already seen in the spectral density that there is no overlap between the 
state of singly occupied wells and a state of a doubly occupied well. From this point of view we expect 
a similar behavior as found for the Hubbard model. However, there is the additional feature that we 
can tune continuously the local atom-atom coupling 7 from an attractive to a repulsive interaction. In 
this way we also reach a degeneracy point at which the effective interaction vanishes (i.e. 7 = 0). The 
existence of only three peaks in the spectral density of Fig. 5 explains the fact that the spin imbalance 
is characterized by only two frequencies, i.e., the difference between the dominant energy level and the 
other levels. 

For <7 = we get the Fermi-Hubbard model without phonon excitations which corresponds with the 
above mentioned experiment. In this case, if the initial state is | '\,\), the dynamics of spin imbalance 



is characterized by two frequencies, U /2 ly/ (AJ/U) 2 + 1 ± lj (sec [25, 27]). The corresponding spin 

imbalance for U — 2 is plotted in Fig. 6. 

For nonzero coupling g the dynamics is affected by the presence of phonon excitations but it is still 
characterized by two frequencies, which are the differences between the second and the first and the second 
and the fourth curves in Fig. 3, respectively, i.e., lo\ = A2 — A3 and uji = A2 — A4. Consequently, at the 
crossing the dynamics of spin imbalance shows only one frequency since the differences between the curves 
in this case are equal, as depicted in Fig. 7. For larger couplings g two frequencies component appear again 
as it is shown in the right panel of Fig. 6. Thus measuring the difference between the two frequencies, 
A = |cj2 — ^i| provides a method to detect the crossing point experimentally: A vanishes at the crossing 
point as shown in Fig. 9 and the spin imbalance is characterized by only one frequency. Increasing the 
coupling g further leaves a low frequency component with almost full amplitude and additional high- 
frequency modulation with small amplitude. This is depicted in Fig. 8. The low frequency component 
oscillates with the frequency u-2 ~ 4r 2 /|7|, where r is given in Eq. (6) and 7 = U — 2g 2 /ujq. Thus LJ2 — > 
as g — > 00, which is a direct consequence of the polaron effect. 




(23) 
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FIG. 7: Spin imbalance Ni^ft) at the level crossing g = 7.06541 for J = 1, U = 2 and ojo = 50. In this case only 
one mode with frequency wi ~ 1.96 contributes to the dynamics. 



VI. DISCUSSION AND CONCLUSION 

At low energies, the restricted model with at most one phonon excitation per well has qualitatively 
the same behavior as the model with many phonon excitations, described by the effective Hamiltonian in 
Eq. (5). This is presented in Figs. 2 and 3, where the four lowest levels are plotted for both cases. The 
main difference, however, is that the Hilbert space of the model with many phonon excitations is much 
larger. Consequently, there are many excited states with energies higher than those shown in Figs. 2, 3. 
However, these states are not considered here because of their high energies. Due to the matrix elements 
K\ and K2 of the effective Hamiltonian in Eq. (11), these higher levels are closely related to harmonic 
oscillator levels with frequency loq. 

Without tunneling (i.e. J = 0) there is a change of the ground state from single occupancy of the wells 
(weak coupling g) to double occupancy of one well (strong coupling g). This reflects the sign change of 
the effective coupling 7 = U — 2g 2 /luq. In the presence of tunneling (i.e. J ^ 0) the ground state, given 
by the coefficients of Eq. (15), changes smoothly upon a change of the coupling g. Its energy is the lowest 
solution of Aa(z) = z, where A3 (2;) is defined in Eq. (14). There is a transition due to the crossing of the 
first and second excited level though (cf. Fig. 3). 

In conclusion, we have studied an atomic mixture of two heavy atoms and two light spin-1/2 fermionic 
atoms in a double- well potential, where the heavy atoms are subject to local harmonic oscillator potentials 
of the wells. This is modeled using the Holstein-Hubbard Hamiltonian, which is the simplest system that 
mimics the presence of phonons in a solid. We have applied the recursive projection method, which 
reduces the complexity of the full Hilbert space and leads to an effective fermionic Hamiltonian. We have 
found a transition for the the light fermions from singly occupied wells to doubly occupied wells as the 
coupling between heavy and light species is increased. This transition is manifested by the crossing of the 
second and third eigenvalue of the effective Hamiltonian. Moreover, the coupling between the light and 
the heavy atoms renormalizes the tunneling of light fermions between wells, which reflects the polaron 
effect. The dynamics is dominated by a spectral density with three peaks. This implies for the spin 
imbalance dynamics of the light atoms a periodic behavior with two characteristic frequencies. These 
frequencies coincide at the crossover of the two lowest excited states. Thus the oscillating behavior of the 
spin imbalance can be used to detect the crossing point experimentally. 

Acknowledgments 

This work was supported by Coordenacao de Apcrfeicoamento dc Pcssoal de Nfvel Superior (CAPES) 
and by the Deutscher Akadcmischer Austausch Dienst (DA AD). 



11 




FIG. 8: Spin imbalance JVi,2(t) for g — 20, J = 1, U — 2 and cjo = 50. Large coupling g leaves a slow 
component with almost full amplitude and an additional high-frequency modulation with small amplitude. The 
slow component oscillates with the frequency UJ2 ~ 4r 2 /|7j. This is the direct consequence of the polaron effect, 

_ 2 / 2 

since it causes the single fermion tunneling decay exponentially, r = Je 9 ' °. 




FIG. 9: Difference A = |cj2 — Wi between the frequencies of the modes contributing to the dynamics of spin 
imbalance. The difference disappears at the crossing point. 
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